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Total reaction cross sections from ^^^Pr(a,a)^^^Pr elastic scattering and a-induced 
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Elastic scattering data for ^*^Pr(Qf,a)^'*^Pr have been analyzed to derive a new energy-dependent 
local potential for the ^*^Pr-a system. This potential is used successfully to predict the cross 
section of the ^''^Pr(Q!,n)^*'*Pm reaction at low energies where new experimental data have become 
available very recently. Contrary to various global potentials, this new potential is able to reproduce 
simultaneously elastic scattering data around and above the Coulomb barrier and reaction data 
below the Coulomb barrier for the ^*^Pr-a system. Reasons for the partial failure of the global 
potentials are explained by intrinsic properties of the scattering matrix and their variation with 
energy. The new local potential may become the basis for the construction of a new global a- 
nucleus potential. 
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I. INTRODUCTION 

The total reaction cross section (Trcac is related to the 
complex scattering matrix Sl = Vl cxp (2i5i) by the 
well-known relation 
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where k — \/2flE^~^Jh is the wave number. i?c.m. is the 
energy in the center-of-mass (cm.) system, and rj^ and 
6l are the real reflexion coefficients and scattering phase 
shifts. (Tl is the contribution of the L-th partial wave to 
the total reaction cross section fTrcac- 

Usually, experimental elastic scattering angular distri- 
butions are analyzed using a complex optical potential. 
At first view, it seems to be a simple task firstly to de- 
termine (Trcac from the analysis of the elastic scattering 
angular distribution and secondly to distribute this cross 
section CTrcac among the open channels (e.g. using the 
statistical model) to predict cross sections of a-induced 
reactions. However, in practice several problems appear. 
There is no unambiguous way to determine reflexion co- 
effcients ijl, phase shifts 6l, or the optical potential from 
a measured elastic scattering angular distribution, and, 
in addition, in most cases angular distributions arc mea- 
sured at relatively high energies whereas reaction cross 
sections should also be known at low energies below the 
Coulomb barrier (this holds in particular for reaction 
cross sections relevant for nuclear astrophysics). Thus, 
typically an ambiguous optical potential has to be ex- 
trapolated down to low energies; as a consequence, con- 
siderable uncertainties have been noticed for the predic- 
tion of a-induced reaction cross sections at low energies, 
in particular for (a, 7) cajrture reactions for targets with 
masses above A « 100 [ll- 
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Very recently, Sauerwein et al. Q have studied the 
-'^*^Pr(a;n)^^^Pm reaction at energies between 10 and 
15MeV, i.e. below the Coulomb barrier. It is shown in 
Q that the calculated ^'^^Pr(a,n)^^*Pm cross section de- 
pends almost exclusively on the a transmission and is 
thus well-suited to test global a-nucleus optical poten- 
tials. It is found that the new experimental data can- 
not be reproduced by any of the widely used global a- 
nucleus potentials which are the energy-independent 4- 
parameter McFadden/Satchler (MCF) potential 0, the 
latest version of the Avrigeanu (AVR) potential with 
many partly energy-dependent parameters [l0| . and the 
energy-independent 6-parameter potential by Frohlich 
and Rauscher (FRR) which is optimized for low-energy 
reaction cross sections [ll|. However, an excellent de- 
scription of the new ^^^Pr(a,n)^^^Pm data is achieved in 
[1] using an energy-dependent modification of the MCF 
potential where a new energy dependence of the depth of 
the imaginary Woods-Saxon volume potential was intro- 
duced to reproduce the new reaction data; the potential 
by Sauerwein et al. is referred to as Sauerwein/Rauscher 
(short "SAR") in the following. 

Contrary to the study in Q which is restricted 
to the analysis of the ^■*^Pr(a,n)^^'*Pm reaction in a 
narrow energy window, the present study considers 
^*^Pr(a,a)^'*^Pr elastic scattering in a wide energy range 
from 19 to 45 MeV. From the fits to the elastic scattering 
angular distributions a new energy-dependent potential 
is derived, and the total reaction cross section (Jreac is cal- 
culated from this potential. (Jrcac is then compared to all 
in the EXFOR data base [l2| available a-induced reac- 
tion data on ^*^Pr. The aim of the present study is thus 
a consistent description of all available elastic scattering 
and reaction data over a broad energy range. 

Elastic ^'*^Pr(a,a)"'^^^Pr scattering data are available 
in literature at 45 MeV [l^. However, these data cover 
only a limited angular range, and they have to be digi- 
tized from Fig. 3 in [l^. The latter leads to uncertain- 
ties which are difficult to estimate. Four ^^-'^Pr(a,a)^^^Pr 
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angular distributions at i^iab = 19.0, 23.97, 32.0, and 
37.7 MeV have been measured by [ij]. The experiment 
has been performed at the XTU Tandem of the INFN 
Laboratori Nazionah di Legnaro. Unfortunately, these 
data have never been published; a partial analysis of the 
data was already shown in an earlier publication p^ . 
These data are available in numerical form (including 
uncertainties) and are thus much better suited for the 
determination of an optical potential by a fitting proce- 
dure. However, as the original data give statistical er- 
rors only (with tiny uncertainties at forward angles), a 
5 % uncertainty has been added quadratically to account 
for unknown systematic errors of the data. Note that 
it is difficult to achieve much smaller systematic uncer- 
tainties for a-scattering experiments (see e.g. [IB, fl7l|). 
These data will be made available to the community via 
EXFOR [13. 

Besides the ^^^Pr(a,n)^'*^Pm data of [§|, the reactions 
i4iPr(a,n)^4pm and "iPr(a,2n)i''3pni have been mea- 
sured in using a stacked-foil activation technique. 
The data cover an energy range from about 15 to 45 McV. 
Although this technique leads to considerable uncertain- 
ties at low energies, these data provide further insight 
into the energy dependence of reaction cross sections and 
the distribution of the total reaction cross section CTroac 
among different open channels. Unfortunately, data for 
^^^Pr(a,7)^'*^Pm are not available in [l2|, and because 
of its lower reaction cross section and unfavorable half- 
life and decay branches of ^^^Pm it was not possible 
to measure the ^^^Pr(a,7)^*^Pm cross section simulta- 
neously with the ^'''^Pr(a,n)-'^^'*Pm cross section in the 
recent high-sensitivity activation experiment of [^. 



II. ANALYSIS OF ELASTIC SCATTERING AND 
THE TOTAL REACTION CROSS SECTION cr^cac 

A. Elastic scattering 

The analysis of the ^''^Pr(a,Q!)^^^Pr angular distribu- 
tions follows closely the procedure outlined in earlier 
work [l^, [l3|. The total potential is composed of the 
nuclear potential with a real and an imaginary part and 
the real Coulomb potential. The different potentials arc 
discussed in the following. 

The real part of the nuclear potential is derived from 
the folding model; the folding potential is modified by a 
strength parameter A w 1.2 — 1.4 and a width parameter 
1/7 ~ 1.0. (Large deviations from w w 1 would indicate a 
failure of the folding model.) To avoid discrete uncertain- 
tics from the so-called "family problem" real potentials 
with volume integrals of about Jn « 320 — 350 MeV fm'^ 
have been selected [l^. (Note that the negative signs of 
volume integrals Jr and J/ are - as usual - omitted in 
the discussion.) 

The imaginary part of the nuclear potential is taken 
in the usual Woods-Saxon paramctrization. It is well- 
known that scattering data at low energies are best de- 



scribed using an imaginary potential of Woods-Saxon 
surface type (see e.g. [13 )• This has been confirmed re- 
cently in a microscopic calculation of the a-nucleus po- 
tential [2^. The same behavior is found in the present 
study where the angular distributions at i?iab = 19.0, 
23.97, and 32.0 McV can be very well reproduced with 
a pure surface imaginary potential. An additional vol- 
ume Woods-Saxon potential in the imaginary part does 
not improve the description of the angular distributions 
at low energies. However, for an excellent description of 
the 37.7 and 45.0 MeV data a combination of volume and 
surface Woods-Saxon potentials for the imaginary part is 
required; at these energies the volume part is even dom- 
inating the imaginary potential. 

The Coulomb potential is taken in the usual 
paramctrization of a homogenously charged sphere with 
a Coulomb radius Rq identical to the root-mean-square 
(rms) radius of the folding potential without width mod- 
ification {w = 1). This avoids uncertainties from the 
otherwise somewhat arbitrary choice of the Coulomb ra- 

1 /3 

dius Rc (often taken as Rc = 1.2, 1.3, or 1.4 fm x AJ ) 
which are non- negligible at least at very low energies [21| . 

The resulting parameters of the potential and the total 
reaction cross section (Troac arc listed in Table IH The fits 
are compared to the experimental angular distributions 
in Fig. [fl 

In addition to the local potential analysis, a model- 
independent phase shift analysis (PSA) has been per- 
formed using the technique of [24]. These phase shift 
fits show the tendency to relatively high and oscillating 
cross sections at backward angles where no experimental 
data are available. Nevertheless, the derived total reac- 
tion cross sections CTroac are close to the results of the 
local potential fit. From the variation of (Troac with the 
fitting parameters (in particular the maximum fitted an- 
gular momentum Lmax) and from the comparison with 
the local potential fit, the uncertainty of (Troac can be es- 
timated to be smaller than 3% in all cases except the 
lowest energy where I estimate an uncertainty of about 
7 %. This is also consistent with a recent analysis of total 
reaction cross sections in [sij . 

As discussed above, the extraction of the total reaction 
cross section (Trcac requires theoretical considerations and 
is thus somewhat model-dependent. Nevertheless, be- 
cause of the small sensitivity of iTreac to the chosen model, 
the total reaction cross section areac can be considered as 
a quasi-experimental quantity. This holds in particular 
for cases where the elastic scattering angular distribu- 
tions cover the full angular range. However, as a word 
of caution, it should be kept in mind that discrepancies 
have been noticed between (Troac determined from elastic 
scattering angular distributions and from (^-transmission 
experiments psi . [23 | . These discrepancies have not been 
fully understood up to now (25| . 

Furthermore, Fig. [1] shows the results of several global 
q;- nucleus optical potentials [9l-[ll| . It is obvious that the 
global potentials cannot achieve the same quality as the 
local potential fit or the phase shift fit. This is not a sur- 
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TABLE I: Parameters of the optical potential and the total reaction cross section (Jreac derived from ^'*^Pr(a,a)^*^Pr angular 
distributions. 



-Elab -Ec.m. 


A 


w 


Jr 


^/?,,rms 


Wv 


Rv 


av 


Ws Rs as 


Ji 






a 

rcac 


Ref. 


(MeV) (MeV) 


(-) 


(-) 


(MeVfm^*) 


(fm) 


(MeV) 


(fm) 


(fm) 


(MeV) (fm) (fm) 


(MeVfm^*) 


(fm) 


(-) 


(mb) 


Exp. 


19.0 18.48 


1.354 


0.994 


345.3 


5.511 








17.1 1.358 0.581 


45.1 


7.451 


2.3 


465±31 




23.97 23.31 


1.270 


1.010 


338.6 


5.604 








27.5 1.410 0.508 


68.2 


7.623 


0.8 


1018±31 


M 


32.0 31.12 


1.452 


0.958 


326.7 


5.314 








36.1 1.290 0.500 


73.9 


7.012 


2.4 


1306±39 


[14J 


37.7 36.66 


1.328 


0.992 


329.3 


5.504 


-17.8 


1.297 


0.157 


3.9 1.586 0.457 


51.7 


6.075 


5.9 


1654±50 


[14J 


45.0 43.76 


1.306 


1.012 


340.0 


5.614 


-18.4 


1.303 


0.107 


6.1 1.567 0.436 


58.5 


6.253 


0.4 


1874±56 


[13] 


low energies 


1.354 


0.994 


345.3 


5.511 








Eq. dU 1.353 0.530 


Eq. @ 


7.361 








a decay 


1.151 


l.OOO' 


308.8 


5.545 





















"from the local potential fit using Eq. JTJ; uncertainties estimated 
from the model-independent phase shift analysis 
'fixed 



prise because the parameters of the global potentials are 
not readjusted to the experimental angular distributions. 
Nevertheless the AVR potential reproduces the angular 
distributions very well. The calculations using the MCF 
potential are also close to the experimental angular dis- 
tributions; however, there is a systematic ovcrcstimation 
of the scattering cross sections at backward angles at low 
energies. A strong overestimation at backward angles is 
found for the FRR potential at all energies; this corre- 
sponds to a significant underestimation of the total re- 
action cross sections iTioac of the FRR potential. It has 
to be pointed out here that the FRR potential was never 
intended to reproduce elastic scattering data above the 
Coulomb barrier. 

The present study does not show results of the po- 
tentials of Kumar et al. [2^ because this potential has 
been optimized for a wide energy range above the new 
i4ipj,j-Q, j^-ji44pj-|^ data; in addition, it has been shown 
in [l^ that this potential cannot reproduce ^^Y{a,a)^^Y 
elastic scattering data at low energies. Unfortunately, 
the latest version s l27l . [28| of the global potential by 
Demetriou et al. [29[ are only published in conference 
proceedings and cannot be used without the authors of 
p7l[28j : I do not intend to show results from the early and 
perhaps out-dated potentials in [2§| . Not yet included arc 
also the predictions from a new regional potential which 
was derived from recent scattering data of the nuclear as- 
trophysics group at Notre Dame university; a publication 
is in preparation [30| . 
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The reduced energy -Ercd takes into account the differ- 
ent heights of the Coulomb barrier in the systems under 
consideration, whereas the reduced reaction cross sec- 
tion (Trod scales the measured total reaction cross section 
CTrcac accordiug to the geometrical size of the projectile- 
plus-target system. A smooth behavior for all Cicd of 
a-induced reactions is found, including the new data for 
^^^Pr-a. The obtained values for (Trod are smaller for 
tightly bound projectiles (a, ^^O) compared to weakly 
bound projectiles (^'^'®Li) and halo projectiles (^He). 
The data are shown in Fig. [5] which is an update of similar 
figures in [3l|, H^j- Up to now, no complete theoretical 
analysis has been presented for the reduced cross sec- 
tions tTred m Fig. [21 and for better readability the data 
points have been connected by lines "to guide the eye" 
(dotted lines in Fig. [2 taken from [s^])- In addition to 
these lines, I show here the result from the local poten- 
tial for ^'*^Pr(Q;,Q!)^'*^Pr with its energy dependence as 
discussed in the next section (Sect. Ill Cp . This calcula- 
tion reproduces practically all data points for a-induccd 
reactions and also data for ^^0-^'^^Ba which is another 
combination of doubly-magic projectile and semi-magic 
target nucleus. 



B. Reduced cross sections CTrcd 

For a comparison of total reaction cross sections for 
different nuclei at various energies, it has been suggested 
to present the data as so-called "reduced cross sections" 
versus "reduced energies" as defined by 

Zjp/jT 



The smooth behavior of all reduced cross sections a^^A 
in Fig.[2]encourages to search for a global potential which 
is able to reproduce drcd and thus the energy dependence 
ot the total reaction cross sections (Treac for a-induced re- 
actions. The present study was restricted to ^"^^Pr but 
an extension to a wider target range is planned for the 
near future. The following procedure is applied to derive 
an energy-dependent potential from the local po- 

tential fits (sec also Table HI . 
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FIG. 1: (Color online) Rutherford normalized elastic scat- 
tering cross sections of the ^'*^Pr(a,a)^''^Pr reaction versus 
the angle in center-of-mass frame. The lines are calculated 
from a local potential fit which is adjusted to the scattering 
data (full red line), from a phase shift analysis (dash-dotted 
lightbrown line) [22|, and from different global a-nucleus po- 
tentials [9l-[Tll|. The experimental data have been taken from 
[l^. [T^. The given energies are iJc.m. in the center-of-mass 
system. 



C. Energy dependence of the potential 

Here I discuss the extraction of a local energy- 
dependent potential for ^^^Pr-a. An estimate of the un- 
certainties for the resulting total reaction cross section 
CTrcac wiU be given later in Sect. IIIIBI 

The energy dependence of the real part is weak. The 
volume integrals Jj^ in Table U increase slightly with de- 
creasing energy; however, at even lower energies the op- 
posite behavior is suggested from dispersion relations. 
Thus, the real part of the potential is simply taken from 
the lowest angular distribution at 19MeV where the 
width parameter w is close to 1.0 and close to the av- 
erage of the other energies. Note that the parameters of 
the real part vary only weakly with energy. This holds 
also for the potential derived from the analysis of the 
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FIG. 2: (Color online) Reduced reaction cross sections 
(Jrod versus reduced energy iJrcd for tightly bound Q-particles 
and ^^O, weakly bound ®'^'*Li projectiles, and exotic ®He. 
(Update of Fig. 4 from [s^l with additional data from 
[3,[i3,llII,lll,H)- The error bars of the new data for "^Pr 
(huge red symbols) are omitted because they are smaller than 
the point size. The dotted lines are to guide the eye. The full 
red line is calculated from the energy-dependent local poten- 
tial for "iPr-a (see Sect.HTCl). 



weak a-decay branch of ^*^Pm (see Sec. Ill D| ). 

The situation for the imaginary part is more difficult 
because the parameters vary with energy, and different 
parametrizations had to be used for the five angular dis- 
tributions under study. For the extrapolation to low en- 
ergies I use a surface Woods-Saxon potential (parameters 
"low energies" in Tablc|l| with the average geometry from 
the angular distributions at 19.0, 23.97, and 32.0 MeV, 
i.e. Rs = 1.353 fm and as = 0.530 fm. Only the three 
lowest energies are considered here because these three 
angular distributions could be described with the same 
paramctrization (surface Woods-Saxon) of the imaginary 
part. Note that the radius parameters Rs and the dif- 
fuscness parameters as do not vary significantly with en- 
ergy, i.e. the shape of the imaginary potential is well de- 
fined from the experimental angular distributions. (A 
significant energy dependence of the shape of the imagi- 
nary part is found only for halo-like projectiles like ^He, 
see e.g. [ssj.) 

The strength Ws of the imaginary potential is derived 
by fitting the imaginary volume integral Jj of the three 
lowest energies using a Fermi-type function (similar to 
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1 + cxp [(£'0 - Ec.m./aE)] 



(4) 



with the parameters J/.o = 74.16MeVfni'^, Eq = 
17.41 MeV, and as = 2.42 MeV; a similar formula like 
Eq. (m holds for the depth Ws of the surface imaginary 
potential with Ws,o = 31.1 MeV. A similar Fermi-type 
function was also used in for the depth of the imagi- 
nary volume Woods-Saxon potential. It is interesting to 
note that the parameters Eq and in Eq. Q are close 
to the values obtained from the adjustment to the new 
"iPr(a,n)"4pm data in Q. 

The widely used Brown-Rho parametrization of the 
imaginary part [35| is given by 



Ji = Ji.o 



{Ec.m. — Eo)'^ 
{Ec.^n. - Eoy + A2 



(5) 



for energies i?c.m. > Eq above the opening of the first 
non-elastic channel at energy Eq and J/ = below Eq. 
Jifi is again the saturation value, and the parameter 
A describes the slope of J/ from zero to its saturation 
value Jifi. Eq is given by the excitation energy of the 
lowest excited state in ^^^Pr which can be pop ulated 
by inelastic scattering: Eq = 0.145McV [H,!!!. The 
Brown-Rho parametrization in Eq. (O is not able to 
reproduce the energy dependence of the imaginary vol- 
ume integrals and thus cannot be used for the extrapo- 
lation of the potential to lower energies. This may - at 
least partly - explain the problems with the prediction of 
the i44Sm(a,7)i48Qd ^^^^^ section Q using the potential 
derived from ^■*^Sm(a,a)^'^^Sm elastic scattering [111 in 
combination with a Brown-Rho parametrization of the 
imaginary volume integral. 

The failure of the Brown-Rho parametrization is shown 
in Fig. [3l Whereas the three parameters of the Fermi- 
type function allow a perfect reproduction of the J/ 
values from elastic scattering, it is impossible to re- 
produce the steep rise of the J/ data between about 
15 and 25 MeV with a Brown-Rho function and fixed 
Eq ~ 0.145 MeV. If the energy Eq of the opening of the 
first non-clastic channel is taken as an additional free 
parameter, the J/ values from elastic scattering can be 
reproduced with J/ = 76.2MeVfm3, Eq = 15.1 MeV, 
and A = 2.8 MeV. However, this result corresponds to a 
vanishing imaginary part already in the energy range of 
the recent ^^^Pr(Q!,n)^*^Pm experiment [1] and thus pre- 
dicts that the total reaction cross section vanishes. It is 
clear that the Brown-Rho parametrization of the imagi- 
nary volume integral J/ cannot be used for the prediction 
of a-induced reaction cross sections of ^"^^Pr. 



D. a-decay of ^*^Pm 

a-decay provides a further opportunity to test a- 
nucleus potentials at low energies. For ^^^Pm (as in most 
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FIG. 3: (Color online) Energy dependence of the imagi- 
nary volume integral Ji using the Fermi-type parametriza- 
tion of Eq. Q (full red line) and Brown-Rho parametriza- 
tions of Eq. (O with fixed Eo = 0.145 MeV (dotted blue 
line) and Eo = 15.1 MeV (dashed green line) with Eo ad- 
justed to the scattering data. The energy range of the recent 
^'*^Pr(Q,n)^''*Pm experiment @] is indicated by vertical lines. 
Further discussion see text. 



other cases) a-decay is dominated by the AL = tran- 
sition from the ^''^Pm ground state to the i^ipj. ground 
state; thus, analysis of a-decay properties provides infor- 
mation on the potential for the L = partial wave (see 
also Sect. |TV] for the energy-dependent relevance of dif- 
ferent partial waves for Creac )- The mass range around 
A 150 has been studied in [3a| using the same type of 
folding potentials as in this work. Although the a-decay 
half-lives vary over many orders of magnitude between 
"SQd, i46Sm, and ^''''Nd (with their N = 82-daughters 
"4Sm, i42Nd, and "°Ce), it has been found in [39[ that 
the preformation of the a-particle in the decaying nucleus 
is between about 10 and 20 % within this model. 

i45pj^-^ Tfj\x\i its N ~ 82-daughter i^^Pr has a weak a- 
decay branch of (2.8 ± 0.6) x lO'^ @, |4l|which to- 
gether with the half-life of T1/2 = 17.7 ± 0.4 y [3i,|43,|42| 
translates to a partial a-decay half-life T^j2 = (2.0 ± 
0.4) X lO^^s. The a-decay Q-value is Qa = 2322.2 ± 
2.6 keV [11]. To repeat the a-decay calculations of [11] 
for ^■^^Pm, a real folding potential has been calculated at 
extremely low energies (labelled "a decay" in Table H]) , 
and the a-decay half-life has been calculated using the 
semi-classical model of [11] . From the ratio between the 
calculated half-life T^jl''^'' = 2.49 x lO^^s and the experi- 
mental partial half-life a preformation of P = 12.5±2.7% 
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is determined which is within the range of 10 — 20 % for 
the neighboring a emitters with N = 82 daughters. Al- 
though based on almost 50 years old data for the half-life 
[4^ and the a branching [4l|, this result nicely confirms 
the close relationship between the various = 82 nuclei 
including ^*'*Sm. 



III. a-INDUCED REACTIONS ON "^PR 

As the optical potential for ^^^Pr-a is completely fixed 
from the above procedure in Sec. Ill CI the calculation 
of the total reaction cross section (Ticac is straightfor- 
ward and does not require any further parameter adjust- 
ment to experimental reaction data. First, the obtained 
o'rcac(^') is Converted to the reduced cross section (T,.od 
and compared to the (Tred data for various projectiles and 
targets (full line in Fig. [2|). It is obvious that the gen- 
eral behavior of the erred vs. E'red energy dependence is 
very nicely reproduced at least down to E^ed ~ 0.8 Me V 
corresponding to i?c.m. ~ 14MeV for the ^^^Pr-a system 
under study in the work. 

Next, the result for (Troac is shown in Fig. |4] in a 
wide energy range, and it is converted to the astro- 
physical S-factor in Fig. [5] for comparison with the new 
^'*^Pr(Q;,n)^^'*Pm data of j^] for energies between 10 
and 15MeV. The new experimental S-factor data de- 
crease with energy by about a factor 5 from S{E) « 
3 X 10^6 MeVb below 11 MeV to S{E) « 0.6 x 10^6 MeVb 
at 14.5 MeV. An excellent reproduction of the new data 
of [1| is found. 



A. Total reaction cross section a^cd and cross 
section of the ^''^Pr(a,n)^^''Pm reaction 

The total reaction cross section (Troac is given by the 
sum over all non-clastic channels, i.e. it includes inelastic 
scattering, fusion, and all transfer channels. This (Trcac 
has to be distributed among the open channels (thresh- 
olds indicated in Fig. |4|) using e.g. the statistical model. 
For the particular case of the ^'*^Pr(Q;,n)^^^Pm reaction in 
the energy range of the recent experiment Q as shown 
in Fig. m it is found that (Treac is almost identical to 
the {a,ii) cross section because of the dominating neu- 
tron emission channel (except very close above the (Q;,n) 
threshold at 10.2 MeV). The proton emission from the 
compound nucleus ^^^Pm is strongly suppressed by the 
Coulomb barrier, and the ratio of (Q;,n) over {a,j) cross 
sections is large in this energy window. Using the stan- 
dard NON-SMOKER parameters [H] (a-potential from 
, nucleon potentials from |45|. 7-ray strength function 
from [4^, level density from |47|) the ratio is sa 4 at the 
lowest energy of Q and exceeds 20 at 11.5 MeV; i.e. ne- 
glecting the weak other open channels like (ajan), the 
total reaction cross section (Troac has to be reduced by 
less than 5 % to obtain the {a,n) cross section above 
11.5 MeV and by w 25% for the two lowest data points. 



The standard TALYS [Hi calculation (of-potential [11], 
nucleon potentials [soj , 7-ray strength function [5lj , level 
density [H, [sll) predicts a ratio of « 10 between {a,n) 
and (q;,7) at the lowest energy of Q and more than 100 
at energies above « 12 MeV; i.e., the {a,n) cross sec- 
tion does not deviate by more than 1 % from (Troac above 
12 MeV, and the deviations at the lowest energies of Q 
never exceed 10%. Thus, because only minor differences 
are predicted between the (a,n) cross section and (Troac 
in the energy range of Q in all calculations [H, \^ , I re- 
strict myself to the presentation of (Treac in Figs. H] and [5l 
This avoids any uncertainties from other sources (mainly 
from the neutron potential and the 7-ray strength func- 
tion) which may be present at the lowest data points of 
Q. The minor differences between the ^^^Pr(Q;,n)^^'*Pm 
cross section and (Troac can be seen from a comparison of 
the calculations in Fig. [5] of this work and Figs. 8 and 9 
of Q (the same scale has been chosen in all these figures 
for simple comparison). 

From Fig. 2] it can be read that (Troac is well repro- 
duced from all potentials under study at higher energies 
above 25 MeV (except the FRR potential which under- 
estimates areac at all energies). Discrepancies become 
visible at lower energies, see Fig. [5j The MCF poten- 
tial shows an incorrect energy dependence and strongly 
overestimates the data at lowest energies. The AVR po- 
tential reproduces the experimental data in the higher 
energy range, but underestimates the data at lower en- 
ergies. The energy dependence calculated from the FRR 
potential is also incorrect, and data at low (high) ener- 
gies are over- (under-) estimated. An energy-dependent 
potential was found in 0] which is able to reproduce the 
new ^^^Pr(a,n)^^^Pm data over the whole measured en- 
ergy range. The prediction of the new energy-dependent 
potential from the present study which was adjusted only 
to elastic scattering data is very close to the result of [1] , 
but does not require any adjustment to reaction data. 

It is interesting to note that the energy of the low- 
est data point in Fig. O corresponds to a reduced energy 
Ercd ~ 0.62 MeV which is close to the lower end of Fig. [2l 
The reduced cross sections (Trod at reduced energies be- 
low E'rod = 1 MeV are thus almost identical for i^ipj. 
+ a and ^^O -I- -'^'^^Ba; in both cases the projectile is a 
doubly-magic nucleus and the target is semi-magic with 
N = 82. 



B. Sensitivity of the total reaction cross section a^cd 
to variations of the parameters of the local potential 

It is difficult to provide a precise error bar for the cal- 
culated total reaction cross section (Troac in Figs. [3] and O 
but an estimate of the uncertainties can be given by rea- 
sonable variations of the parameters of the local potential 
which are the geometry (in particular of the imaginary 
part) and the energy dependence. The results of the fol- 
lowing sensitivity study are shown in Fig. [B] The local 
potential calculation (full red line, identical to Fig. [5]) is 
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FIG. 4: (Color online) Total reaction cross sections (Jreac for a-induced reactions on ^*^Pr from elastic scatte ring (red dots) 
and cross sections of the ^*^Pr(a,n)^''''Pm (open and full triangles [1, [l3|) and ^*^Pr(«,2n)^'''^Pm (full squares [isl ]) reactions. 
The data of [3 are connected by thin lines to guide the eye. Calculations with different potentials [3-[Tll] agree well at higher 
energies (except [ll[), but disagree at lower energies (more details see also next Fig. The arrows on top indicate the 
thresholds for various reaction channels. 



the reference for the sensitivity study. 

The real part of the potential is nicely constrained 
from the folding calculation. In addition to the scattering 
data, further information on the real part of the potential 
at very low energies can be extracted from the analysis 
of the a-decay of ^^^Pm (see Sect. Ill D|) . The parameters 
of the potential for the a-decay calculations remain close 
to the real part of the scattering potential. The volume 
integral Ju is about 10 % lower than the value found for 
the 19MeV scattering data. If this lower real part of 
the potential is used for the calculations of cTroac instead 
of the 19MeV real part (as derived in Sec. Ill C|) . (Trcac 
shows a slightly different energy dependence but does 
not change by more than about 40 % in the energy range 
under study (dashed fuchsia line in Fig. ^ . Only at very 
low energies below lOMeV (i.e., below the shown energy 
range of Fig. [B]), the lower real potential with its resulting 
higher effective Coulomb barrier results in a lower Urcac. 
But in any case the deviations remain below 50%, thus 
confirming that the real part of the potential is relatively 
well-defined and does not lead to big uncertainties in the 
calculation of (ircac in the energy range under analysis. 

It should be noted that the a-decay potential is ad- 
justed at the decay energy of 2.3 MeV far below the en- 
ergy range under study; thus, using the a-decay potential 
for the calculation of (Trcac is an extremely careful esti- 
mate for the uncertainty of CTreac on a variation of the 
real potential. I do not show calculations with increased 
strength of the real part because volume integrals sig- 
nificantly above 350 MeV fm'^ have not been observed in 
a-scattcring of semi- magic nuclei [l^ [l^ [l^ [5^ . This 
finding is supported theoretically by the fact that disper- 



sion relations lead to a reduction of real potential at very 
low energies (see e.g. Fig. 11 of [IB] or Fig. 12 of [s^]). 

The influence of the imaginary part on the calculated 
total reaction cross section CTreac is significant. In par- 
ticular, a reasonable energy dependence is essential for 
the reproduction of the experimental data. If the en- 
ergy dependence of the imaginary part is ignored and 
the saturation value J/^o is used instead of the energy 
dependence of Eq. then CTreac is overestimated with 
increasing discrepancy to the experimental data at lower 
energies (green dotted line in Fig. ^ . A similar behavior 
has been found for the energy-independent MCF poten- 
tial. However, small variations of the imaginary strength 
result only in minor modifications of CTreac : a reduction of 
Ji{E) by 10% leads to slightly reduced CTreac with very 
similar energy dependence as the reference calculation 
(dash-dotted aquamarine). 

Next the geometry parameters of the surface Woods- 
Saxon potential in the imaginary part were taken from 
the three angular distributions at 19, 24, and 32 MeV (see 
Table m instead of their average values Rs = 1.353 fm 
and as = 0.530 fm. The depth Ws has been adjusted to 
the same volume integral in Eq. ^ as in the reference 
calculation. The geometry from the 19 MeV angular dis- 
tribution leads to somewhat larger CTroac which is a con- 
sequence of the larger 05 = 0.581 fm (short-dashed blue). 
The result from the 24 MeV geometry is almost identi- 
cal to the reference calculation (dash-dotted magenta), 
and the result from the 32 MeV geometry is somwehat 
smaller than the reference (dash-dotted golden). The 
uncertainty of the geometry of the imaginary potential 
leads to uncertainties of CTreac of the order of 30%. It 
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FIG. 5: (Color online) Same as Fig. 3] but shown as astro- 
physical S-factor in a narrow energy window around the new 
experimental data of 0| . For simpUcity, the averaged S-factor 
at the mean energy (see Table III of [gl) is shown here. The 
local potential (full red line) and the SAR potential (dash- 
dotted golden) reproduce the data well, but the global poten- 
tials fail to reproduce the energy dependence (MCF: dotted 
green; AYR: dashed blue; FRR: dash-dotted dark- magenta). 

has to be pointed out that such a small uncertainty can 
only be achieved because the geometry parameters are 
well-defined from several angular distributions at ener- 
gies around and slightly above the Coulomb barrier. 

Summarizing the above sensitivity study, it is found 
that the uncertainty of (Troac from the present local poten- 
tial is much smaller than the variations between the dif- 
ferent predictions from global potentials. The influence 
of the geometry of the imaginary part is not dramatic 
as long as the geometry is well-defined from low-energy 
scattering data. However, the energy dependence of the 
strength J/ of the imaginary part is an essential ingredi- 
ent for the prediction of reaction cross sections below the 
Coulomb barrier. 



IV. DISCUSSION 

Let me start the discussion with a few general remarks 
on statistical model calculations. In the statistical model 
the cross section for a reaction is given by the product 
of the compound formation cross section and the decay 
branching of the decaying compound nucleus into the 
particular channel. The full formalism can be found e.g. 
in [211, 13] ■ The formation cross section is calculated from 



FIG. 6: (Color online) Sensitivity of the astrophysical S- 
factor to variations of the potential parameters. The full red 
line is identical to the previous Fig. [S] A reduced real po- 
tential (derived from a-decay, see Sect. ITT^ shows a slighly 
different energy dependence (dashed magenta). Neglecting 
the energy dependence of the imaginary part leads to strong 
overestimation of the experimental cross sections (green dot- 
ted), whereas small modifications of the imaginary part (10% 
reduction of the strength or variation of geometry parameters) 
do not affect the calculated Ureac strongly. Further discussion 
see text. 



transmission coefficients using the optical potentials, i.e. 
it is the total reaction cross section cTrcac in Eq. ([T]). The 
decay branching is also calculated from transmission coef- 
ficients for particle channels and from the 7-ray strength 
function for the photon (capture) channel. It is obvi- 
ous that the correct reproduction of the total reaction 
cross section (Trcac is the basic prerequisite for a success- 
ful prediction of any reaction cross section in the statis- 
tical model and should thus always be discussed first in 
presentations of statistical model calculations. Unfortu- 
nately, this is not always the case. 

The total reaction cross section drcac is then dis- 
tributed among all open channels. It is important to 
study the influence of all open channels and their rel- 
evance for the decay branching of the compound nu- 
cleus. This is nicely illustrated in Fig. 21 Close above 
the opening of the (a,2n) channel around 17 MeV, the 
(a,2n) channel becomes stronger than the (a,n) chan- 
nel which was dominant from close above its threshold 

at 10.2 MeV up to 17 MeV. If an open channel, e.g. the 
"iPr(Q;,2n)i43pj^ 

channel above 17 MeV, is not taken 
into account in a statistical model calculation, then all 
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the calculated ^'^^Pr{a,X) cross sections above 17MeV 
must be overestimated. However, such a shortcoming of 
neglected open channels may be partly compensated by 
using special potentials for the particular {a,X) reaction 
under study. 

After these general remarks let me be more specific 
for a-induced reactions on ^^^Pr. The total reaction 
cross section cTioac is well understood and has relatively 
small uncertainties for energies above the Coulomb bar- 
rier. This can be seen from the comparison of predicted 
Croac from the various global potentials [9l-IT]]| to the ex- 
perimental result from the analysis of the elastic scatter- 
ing angular distributions (see Table HI)) . All potentials 
(except the FRR potential) reproduce the experimental 
freac within minor uncertainties of typically a few per 
cent. 

TABLE II: Experimental total reaction cross sections (Trcac (in 
mb) derived from ^*^Pr(Q,Q)^^^Pr elastic scattering angular 
distributions (from Table|l]), compared to the phase shift anal- 
ysis ("PSA"), predictions from the global a-nucleus poten- 
tials MCF, AVR, and FRR and the energy-dependent 
potential SAR of d which was adjusted to reproduce the 
^''^Pr(Q,n)^*'*Pm cross section at low energies. 

^c.m. (MeV) exp. PSA MCF AVR FRR SAR 



18.48 


465±31 


433 


425 


437 


197 


424 


23.31 


1018±31 


1003 


957 


993 


648 


957 


31.12 


1306±39 


1310 


1452 


1487 


1125 


1452 


36.66 


1654±50 


1654 


1662 


1696 


1333 


1662 


43.76 


1874±56 


1873 


1841 


1872 


1514 


1841 



The excellent reproduction of (Trcac for most global po- 
tentials at energies above the Coulomb barrier can be 
explained by a look at the reflexion coefficients r]L and 
the contribution ctl of the L-th partial wave to (Xreac (see 
Fig. [7] for the example of the i^c.m. = 23.31 MeV data). 
For small angular momenta L the reflexion coefficients 
are close to rjL ~ 0; thus, ctl increases linearly propor- 
tional to 2L + 1 (dashed line in Fig. [7|). All global poten- 
tials with a reasonable real part and a sufficiently strong 
imaginary part, i.e. the MCF, AVR, and SAR potentials, 
predict the same itl for small L. For large L (correspond- 
ing to large impact parameters) the 77^ approach unity, 
and the vanish (i.e., partial waves with large angular 
momenta are practically not absorbed by any imaginary 
potential with a limited radial range). Differences in 
appear only for a few partial waves between L w 8 and 
L « 14 (for the Ec.m. = 23.31 MeV case) where the dif- 
ferent potentials predict different 77^ and thus different 
aL- It becomes obvious that relatively small differences 
in r]L for a very limited number of partial waves around 
L « 10 cannot lead to major differences for the total reac- 
tion cross section (Trcac which is simply given by the sum 
over all cr^. The good reproduction of the total reaction 
cross section (Treac from different global potentials is thus 
simply a consequence of intrinsic properties of the scat- 
tering matrix as long as the energy exceeds the Coulomb 
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FIG. 7: (Color online) Reflexion coefficients r\L (upper, (a)) 
and contribution ol of the L-th partial wave to the total 
reaction cross section (Treac (lower, (b); linear scale). The thin 
dotted lines connect the data points to guide the eye. Further 
discussion see text. 

The situation changes dramatically at lower energies 
far below the Coulomb barrier. Here only the lowest 
angular momenta (below L w 5) arc affected, show 
< 1, and thus contribute to cr^; but even these 77/, 
remain much larger than zero and do even approach 
77L w 1 for very low energies (see Fig. [8] for energies 
11 MeV < i?c.m. < 15 MeV, i.e. the energy range of the 
new "iPr(a,n)"4Pm data of dl). Here the trivial 2L + 1 
proportionality of the vanishes. Now the relevant cr^, 
depend sensitively on details of the potential. In particu- 
lar, the strength of the imaginary part for large radii has 
strong impact on the resulting cross sections ol- Under 
these conditions it becomes mandatory that the imagi- 
nary potential has the proper geometry (e.g. fixed from 
scattering data at slightly higher energies) and the proper 
strength (using a realistic energy dependence, e.g. from 
Eq. (|4]), which can be adjusted to scattering data and/or 
properly chosen reaction data). As pointed out in SccHIl 
the real part of the potential has only minor uncertain- 
ties in shape and strength because it is calculated from 
a folding procedure; furthermore, the analysis of a-dccay 
properties provides a test for the L = potential (see 
Sect.[IlD|. 

The above discussion explains why all global potentials 
are able to reproduce total reaction cross sections (Trcac 
above the Coulomb barrier, but may fail to reproduce 
(Trcac bclow tlic Coulomb barrier. Thus, it is not surpris- 
ing that the global potentials show a significant scatter in 
the predictions of the new experimental ^'*^Pr((a;,n)^^'*Pm 
data of [1] (see Fig. [3). 

For the energy-independent MCF potential the rea- 
son for the discrepancy is obvious. The missing en- 
ergy dependence of the MCF potential leads to an 
overestimation of the imaginary part of the potential 
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FIG. 8: (Color online) Reflexion coefficients rjL (upper, (a)) 
and contribution ctl of the L-th partial wave to the total 
reaction cross section o"rcac (lower, (b); logarithmic scale) for 
low energies between 11 and 15MeV. Note the different scale 
for TjL in Figs. [7] and |8] The thin dotted lines connect the 
data points to guide the eye. The tjl at 11 and 12MeV are 
very close to unity for all L; thus, the 11 MeV data disappear 
behind the 12 MeV data in the upper diagram. Differences 
between the 11 and 12 MeV data become visible only in the 
presentation of which are proportional to (1 — 7?i). The 
2L + 1 proportionality oi for small L (as seen at higher 
energies in Fig. [7]) vanishes completely at energies below the 
Coulomb barrier. Further discussion see text. 



at low energies and thus to an overestimation of the 
^^^Pr(a,n)^'*'*Pni cross section. As expected, the dis- 
crepancy between the MCF prediction and the experi- 
mental data increases with decreasing energy. A simi- 
lar behavior has been found in other cases, see e.g. the 
i'^''Sm(a,7)"8Gd reaction [l|, the ^^^Sn{arfy^'^Te reac- 



tion j3|J, or the i''<'Cd(a,7)"°Sn reaction ^. 

It is difficult to make such general statements on the 
many-parameter AVR potential. This potential has been 
adjusted to a huge data base of elastic scattering and re- 
action data (ssf , and excellent agreement has been found 
for many reactions especially in the A w 100 range. The 
AVR potential shows the best agreement of the global 
potentials with the elastic scattering angular distribu- 
tions, and it reproduces the ^^^Pr(a,n)^''^Pm data well 
at higher energies above 13 MeV. However, it underesti- 
mates the data at lower energies significantly. 

The FRR potential has been optimized for reac- 
tion data at low energies. Its parameters are energy- 
independent and close to the MCF potential in most 
cases. Similar to the MCF case, because of its energy 
independence it cannot be expected that the FRR poten- 
tial is able to reproduce simultaneously elastic scattering 
data above the Coulomb barrier and reaction data below 
the Coulomb barrier. Despite the big success in predict- 
ing reaction cross sections at low energies, the simple 
FRR potential cannot be the basis for a global potential 
in a broad energy range. 



The MCF potential with an additional energy depen- 
dence of the depth of the imaginary potential has been 
used in to reproduce their new ^^^Pr(a,n)^''^Pm data. 
Although this is in principle the best way to improve the 
MCF potential, there are two disadvantages of the SAR 
potential. Obviously, reaction data are required to fit the 
energy dependence of the SAR potential. And in addi- 
tion, the chosen underlying MCF potential uses a volume 
Woods-Saxon imaginary part which has been found - ac- 
cording to recent studies [13, [H, - to be probably 
not fully adequate for low energies below the Coulomb 
barrier; here surface Woods-Saxon potentials should be 
preferred. Consequently, the SAR potential is not able 
to reproduce the elastic scattering angular distributions. 
It turns out that the angular distributions of the SAR 
and MCF potentials are almost identical above 25 MeV 
because the MCF and SAR potentials are almost iden- 
tical there. Even at 19 MeV where the imaginary depth 
of the SAR potential is about 25 % lower than the MCF 
depth, the calculated angular distributions of MCF and 
SAR agree within a few per cent, i.e. within a linewidth 
in Fig. [U therefore, the SAR potential is not shown in 

Fig.m 

The energy-dependent potential of this work has been 
derived from elastic scattering data. Local optical poten- 
tials have been derived from the available angular distri- 
butions between 19 and 45 MeV by fitting the parameters 
of the real part (strength and width of the folding poten- 
tial) and the imaginary part (Woods-Saxon parametriza- 
tion) . Because of the small variation of the found param- 
eters at the different energies (see Table |l|, these param- 
eters could be combined to derive a common potential 
with an energy-dependent depth of the imaginary surface 
potential. This common potential maintains the good re- 
production of the elastic scattering data and the derived 
total reaction cross sections CTicac- In addition, it turns 
out that this potential is able to predict the cross section 
of the ^'''^Pr(a,n)^''^Pm cross section at lower energies, in 
particular in the energy range of the new experimental 
data of [il between 10 and 15 MeV (see Fig. [5]). Thus, 
the new potential is able to describe ^''^Pr(a,a)^"'"'^Pr 
clastic scattering data at energies around and above the 
Coulomb barrier simultaneously with ^*^Pr(a,n)'^*^Pm 
reaction data below the Coulomb barrier. 

The parametrization of this new potential, i.e. a dou- 
ble folding potential in the real part and a surface 
Woods-Saxon potential with fixed geometry and energy- 
dependent depth, should be tested in further cases in 
forthcoming work. For a detailed study of a particu- 
lar target nucleus, several elastic scattering angular dis- 
tributions are required to fix the shape and the energy 
dependence of the imaginary potential. However, the 
found geometry parameters in the present ^^^Pr-a case 
{Rs = 1.353 fm and as = 0.530 fm) are close to standard 
values; it should be possible to apply these parameters 
to other target nuclei. The main problem will be the de- 
termination of the energy dependence of the depth of the 
imaginary surface potential. Here considerations simi- 
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lar to [1, [13, will be helpful where parameters of the 
imaginary part have been put in relation to the Coulomb 
barrier. It is obvious that the predictions from such a 
potential have to be compared to experimental reaction 
data below the Coulomb barrier which are - despite sig- 
nificant improvement in the recent years - still not widely 
available. 

The previous discussion may also be summarized to 
provide a general recipe to find a widely useful potential. 
First of all, the potential has to reproduce the total re- 
action cross section a^eac at energies above the Coulomb 
barrier. However, because of the 2L + 1 proportionality 
of aL for small L, this requirement is almost trivial and 
can be fulfilled by many potentials. Next, among the 
many a-nucleus potentials which are able to fulfill this 
first condition, these potentials have to be selected which 
are in addition able to reproduce the elastic scattering 
angular distributions. This ensures the correct strength 
and geometry of the chosen potential (in particular of 
the imaginary part) which is essential for the prediction 
of (Treac bclow the Coulomb barrier. This conclusion is 
also a strong motivation to extend the measurements of 
elastic scattering angular distributions. 

V. SUMMARY AND CONCLUSIONS 

A new local potential for the system ^'^^Pr-a has been 
derived from elastic scattering angular distributions. The 
derived potential is able to predict the total reaction cross 
section (Troac which is dominated by the ^'*^Pr(Q;,n)^^'*Pm 
cross section in the energy range of the new experimental 
data of is]. Thus, the new potential for ^^^Pr-a is able for 
the first time to describe simultaneously elastic scattering 
data around and above the Coulomb barrier and reaction 
data below the Coulomb barrier. Such a simultaneous 
description was not achieved and/or aimed for in earlier 
work using the global potentials of [ol-fTlj or the local 
potential suggested in 

Reasons for the partial success at higher energies and 
failure at lower energies of global potentials are carefully 
analyzed by studying the reflexion coefficients rj^ and 
the contribution (Tl of the L-th partial wave to the to- 



tal reaction cross section CTreac- It is found that above 
the Coulomb barrier for small angular momenta are 
predicted correctly by most global potentials because of 
a simple 2L + 1 proportionality which is almost inde- 
pendent of details of the optical potential. This results 
in minor differences for the total reaction cross section 
Crcac which havc to arise from a few partial waves around 
L w 10 for the ^^^Pr-a case under study in this work. 
However, at low energies below the Coulomb barrier this 
simple 2L + 1 proportionality of for small L vanishes, 
the ctl become sensitive to details of the potential, and 
consequently predictions for CTreac from global potentials 
show a huge spread. Precise scattering data at differ- 
ent energies around and close above the Coulomb barrier 
arc required to determine a potential (in particular the 
shape and strength of its imaginary part) which is able to 
predict total reaction cross sections CTroac not only above, 
but also below the Coulomb barrier. 

Because the new local potential has a well-defined ge- 
ometry derived from elastic scattering, its extrapolation 
to low energies and its predictions for the ctl and the 
resulting CTioac should be more reliable than predictions 
from other global potentials. The parametrization of the 
new local potential may finally become the basis for a 
new global potential to solve or at least reduce the long- 
standing problem of a-nucleus potentials at low energies 
below the Coulomb barrier and the resulting uncertain- 
ties for the prediction of a-induced reaction cross sections 
at astrophysically relevant energies. 
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